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Stellar explosions such as novae and supernovae produce most of the heavy elements in the 
Universe. Although the onset of novae from runaway thermonuclear fusion reactions on the 
surface of a white dwarf in a binary star system is understood 1 , the structure, dynamics, and 
mass of the ejecta are not well known. In rare cases, the white dwarf is embedded in the 
wind nebula of a red- giant companion; the explosion products plow through the nebula and 
produce X-ray emission. Early this year, an eruption of the recurrent nova RS Ophiuchi 2 ' 3 
provided the first opportunity to perform comprehensive X-ray observations of such an event 
and diagnose conditions within the ejecta. Here we show that the hard X-ray emission from 
RS Ophiuchi early in the eruption emanates from behind a blast wave, or outward-moving 
shock wave, that expanded freely for less than 2 days and then decelerated due to interaction 
with the nebula. The X-rays faded rapidly, suggesting that the blast wave deviates from 
the standard spherical shell structure 4-6 . The early onset of deceleration indicates that the 
ejected shell had a low mass, the white dwarf has a high mass 7 , and that RS Ophiuchi is a 
progenitor of the type of supernova integral to studies of the expansion of the universe. 

On February 12, 2006, RS Ophiuchi (RS Oph) was discovered in outburst near the historical 
maximum optical brightness of 4.4 magnitudes 8 - more than 4 times brighter than the threshold for 
naked-eye visibility. The Rossi X-Ray Timing Explorer (RXTE) satellite 9 first detected strong hard 
X-rays from RS Oph - the strongest ever observed from a white dwarf - at peak X-ray brightness 
on outburst day 2 with the All Sky Monitor camera. RXTE subsequently measured the X-ray 
spectrum with the Proportional Counter Array (PCA) on outburst days 3, 6, 10, 14, 17, and 21 
(Fig. 1). During this time, the binary was in the phase of its orbit in which the white dwarf and red 
giant appear side-by-side from the perspective of the earth 10 ' 11 . 



During the first three weeks of the outburst, the total X-ray luminosity from the hottest X-ray 
emitting gas fell from ~ 300 L Q to ~ 30 L Q (integrated from 0.5 to 20 keV, and taking a distance 
of 1.6 kpc 1213 ). From day 6 onward, the X-ray emission from the hottest gas decreased as t~ 5 ^ 3 , 
where t is the time since the beginning of the outburst (Fig. 2a). After rising to maximum in about 
2 days, the emission measure (n 2 x V, where n is the density of radiating electrons and V is the 
X-ray emitting volume) remained roughly constant at ~ 10 58 cm~ 3 until day 6, after which time 
it fell as r 4/3 . For comparison, the X-ray luminosity from hot gas in classical novae, where the 
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white dwarf is not embedded in a dense nebula, is always below 30 ' 15 , and the classical nova 
with the most comprehensive X-ray observations 16 (VI 974 Cyg) took more than 100 days to reach 
a maximum X-ray emission measure of a few times 10 56 cm -3 . As the outburst progressed and the 
X-ray luminosity decreased, the bulk of the X-ray emission from RS Oph shifted to lower energies 
(see Fig. 1). This shift corresponds to a drop in the post-shock temperature. From the time of the 
first spectral observation on day 3, this temperature decreased as t" 2 ^ (Fig. 2b). 

If the X-rays are from circumstellar material heated by the blast wave produced in the nova 
explosion, we can associate the characteristic X-ray temperature with the blast wave speed via the 
standard strong-shock relation. Optical emission-line widths directly measure the motion of either 
gas behind the blast-wave shock that has undergone charge exchange 17 or the ejecta, and these line 
widths 18 support the shock speeds we derive. Since v oc T 1/2 , where T is the post-shock plasma 
temperature and v is the blast-wave speed, the rate of temperature decrease indicates that the blast- 
wave speed decreased as t -1 / 3 . This direct measurement of the shock deceleration rate is the first 
for any nova. For an initial shock speed of 3,500 km/s 19 ' 20 , and extrapolating backwards from our 
inferred shock speeds, we conclude that the deceleration began at approximately day 1.7. With the 
measured deceleration, we can obtain the blast wave radius as a function of time. Comparing the 
blast wave radius on days 21 and 27 to radio images of the expanding shock on those days (Rupen, 
Mioduszewski, and J.L.S., in preparation), we confirm the 1.6 kpc distance to RS Oph. 

The standard blast wave theory presupposes that the shock structure is self-similar, or that 
the density and temperature as a function of fractional distance from the white dwarf to the shock 
front are similar to these functions at later times. In the initial ejecta-dominated phase, the ejecta 
expand freely and produce a shock moving at constant velocity. After the expanding blast wave 
has swept up a few times the ejecta mass, it begins to decelerate and enters the Sedov-Taylor (ST) 
phase. After the 1985 outburst of RS Oph, standard supernova blast wave theory was adapted to 
explain the late-time X-ray emission 6 . Models for RS Oph must also account for the fact that the 
nova explosion takes place within a stellar wind whose density far from the binary decreases as 
the inverse of the distance from the binary squared (n oc 1/r 2 , where r is the distance from the 
wind-producing red giant). Whereas the early hard X-ray emission from classical novae is thought 
to arise in the ejecta that has been heated by a reverse shock 21 , the early hard X-ray emission from 
RS Oph comes predominantly from circumstellar material heated by the forward-moving blast 
wave. 

Our observations indicate that the standard self-similar blast-wave model does not explain 
the ejection of material from this nova. Although our discovery that v oc t -1 / 3 agrees with the 
theoretical prediction for a blast wave moving into a stellar wind and the general expansion law 
for a shock wave resulting from a point explosion 4-6 , the observed rate of X-ray fading disagrees 
with predictions. In the current outburst, the initial X-ray flux was almost 100 times brighter 
than expected based on observations and modeling of the late-time X-ray emission in the 1985 
outburst 6 22 . The total X-ray flux for thermal bremsstrahlung emission from a fully ionized gas is 
proportional to T^^rfiV . If the volume of the post-shock emitting region increases as the radius 
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of the blast wave cubed (as in the standard similarity solutions for supernova shells), the X-rays 
should decrease as i -1 . The X-ray emission from RS Oph in the ST phase decayed instead as t~ 5 / 3 . 

The X-ray flux evolution reflects the environment of the binary and suggests how the nova 
explosion deviates from a spherically symmetric, self-similar blast wave. The fast X-ray rise of 
RS Oph compared to classical novae is due to the higher circumstellar density and the faster initial 
shock speed. The observed peak emission measure can be attained in 2 days if the material swept 
up by the blast wave produces thermal bremsstrahlung emission at the immediate post-shock den- 
sity and temperature. The observed X-ray decay rate from day 6 onward is consistent with emission 
from a hot post-shock region whose volume increases roughly as the square of the blast-wave ra- 
dius. Such an emitting-volume evolution could be produced by radiative cooling, which has little 
impact on the early-time shock velocity, but a significant impact on the post-shock temperature 
structure 23 . Alternatively, or in conjunction with radiative cooling, a non-spherical ejecta geom- 
etry such as that inferred from radio observations 24 ' 25 , or a contribution to the X-ray luminosity 
from ejecta heated by a reverse shock, could also steepen the X-ray decline. 

Finally, since the transition to the ST phase occurs when a few times the ejecta mass have 
been swept up, the timing of the onset of deceleration constrains the ejected shell mass indepen- 
dently of all previous estimates. Taking a density inside the binary of 10 9 cm -3 from observational 
constraints on the red-giant mass-loss rate 10 and assuming that the density fall-off does not begin 
until outside the binary, on ~day 2 the shock had swept up only a few times 1CT 7 M of mate- 
rial. Thus, the ejecta mass could not have been much more than 1CT 7 M . Theoretical models of 
nova explosions with recurrence times of 20 years, as in RS Oph, have been constructed for white 
dwarfs with masses of 1.25 and 1.40 M 7 . Whereas the 1.25 M Q models require an ejecta mass 
of ~ 1(T 6 M Q , the 1.40 M Q models eject material with a total mass of 2 x 1(T 7 M . The ejected 
shell mass of ~ 1CT 7 M from RS Oph is thus more consistent with the 1.4 M model. 

The mass accumulation efficiency of white dwarfs in recurrent novae is controversial 7,26 . Our 
conclusion that the white dwarf in RS Oph must be extremely close to the maximum mass shows 
that significant mass accumulation is possible, and that some recurrent novae will explode as Type 
la supernovae. If recurrent novae lead to supernovae la such as the one recently discovered to 
contain hydrogen in its spectrum 27 , sweeping up of the red-giant wind by the nova blast wave could 
explain the cavity surrounding this supernova 28 29 . Our results also have bearing on studies of Type 
II supernovae. The early days of the 2006 outburst of RS Oph show the first two phases of blast 
wave evolution, the analog of which can take hundreds of years for a supernova remnant. They 
also provide the blast-wave deceleration law that is difficult to measure directly in supernovae. Our 
results should motivate modeling that will lead to better understanding of cooling mechanisms as 
well as post-shock and ejecta structure of novae and other stellar explosions. 
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Figure 1 X-ray spectra from the first 3 weeks of the 2006 outburst of RS Oph. These six spectra 
were taken with the PCA instrument on board the RXTE satellite. The abscissa is the energy, E, of 
the detected X-rays. The quantity E x F E is the photon energy times the energy flux density, and it 
provides a measure of the relative amount of energy being emitted across the spectrum. The usable 
energy range for the PCA is approximately 2-25 keV. The spectra were reduced using standard 
data reduction software (the FTOOLS package) and the standard RXTE background models. The 
presence of a blended emission line from H-like and He-like Fe indicates that the X-ray emission 
is optically thin thermal plasma emission; the spectra can all be reasonably well fit (reduced \ 2 
ranging from 0.4 to 1.7 for the 6 observations) with a single-temperature thermal bremsstrahlung 
model plus line emission from Fe and absorption by intervening material. Although RXTE is not 
very sensitive to absorption, which preferentially affects photons with energies less than around 
2 keV, the absorption was high enough on day 3 (Nh = [5.5 ± 1.1] x 10 22 cm -2 , where Nh is 
the column density of neutral hydrogen) to have a significant impact on the low-energy end of the 
spectrum. The absorption lessened in the subsequent observations as the blast wave moved outward 
and the amount of neutral intervening wind decreased. After day 10, we fixed the absorption 
parameter in our fits to the interstellar value 13 ' 22 ' 30 of Nh = 4 x 10 21 cm -2 . Error bars represent 1 
standard deviation. 
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Figure 2 X-ray flux and post-shock plasma temperature as a function of time. We measure the 
time, t, from the start of the outburst, which we take to be February 12.6 U.T. a, X-ray flux in 
the energy range 0.5-20 keV from the hottest plasma, obtained by integrating the best-fit thermal 
bremsstrahlung model, corrected for the effects of absorption by intervening material. This energy 
range was chosen to match typical observing ranges and to include most of the emission from the 
hot plasma. The uncertainty was taken to be 5%, and comes primarily from uncertainty in the 
model parameters, especially interstellar absorption, b, The characteristic X-ray plasma tempera- 
ture determined from fitting the X-ray spectra with thermal bremsstrahlung models, as a function 
of time. Error bars represent 1 standard deviation. The temperature decays as expected for mate- 
rial that has been shock-heated by a decelerating Sedov-Taylor blast wave moving into an oc 1/r 2 
stellar wind. Since RXTE is sensitive to X-rays with energies greater than 2 keV, we obtain the best 
determination of the parameters from the first several observations, when the post-shock plasma is 
hottest. By the sixth observation, the error bars have become too large for these data to be useful 
(see Fig. 1), so the data from the final observation are not included on the plots. 



5 



1. Starrfield, S., Truran, J. W., Sparks, W. M. CNO abundances and hydrodynamic studies of 
the Nova outburst. V - 1.0-solar-mass models with small mass envelopes. Astrophys. J. 226, 
186-202(1978). 

2. Kenyon, S. J. The Symbiotic Stars. (Cambridge: Cambridge University Press, 1986). 

3. Bode, M. R, Ed. RS Ophiuchi (1985) and the Recurrent Nova Phenomenon. (Utrecht: Vnu 
Science Press, 1987). 

4. Sedov, L. Similarity and Dimensional Methods in Mechanics (New York: Academic, 1959). 

5. Chevalier, R. A. Self-Similar Solutions for the Interaction of Stellar Ejecta with an External 
Medium. Astrophys. J., 258, 790-797 (1982). 

6. Bode, M. E, & Kahn, F. D. A Model for the Outburst of Nova RS Ophiuchi in 1985. Mon. 
Not. R. Astr. Soc. 217, 205-215 (1985). 

7. Yaron, O., Prialnik, D., Shara, M. M., & Kovetz, A. An Extended Grid of Nova Models. II. 
The Parameter Space of Nova Outbursts. Astrophys. J. 623, 398-410 (2005). 

8. Narumi, H., et al. RS Ophiuchi. IAU Circ, 8671 (2006). 

9. Jahoda, K., Swank, J. H., Giles, A. B., Stark, M. J., Strohmayer, T., Zhang, W., & Morgan, 
E. H. In-orbit performance and calibration of the Rossi X-ray Timing Explorer (RXTE) Pro- 
portional Counter Array (PCA). Proc. SPIE, EUV, X-Ray, and Gamma-Ray Instrumentation of 
Astronomy VII, O. H. Siegmund, & M. A. Gummin, Eds., 2808, 59-70 (1996). 

10. Dobrzycka, D., & Kenyon, S. J. A new spectroscopic orbit for RS Ophiuchi. Astron. J., 108, 
2259-2266 (1994). 

1 1 . Fekel, F. C, Joyce, R. R., Hinkle, K. H., & Skrutskie, M. F. Infrared spectroscopy of symbiotic 
stars. I. Orbits for well-known S-type system. Astron. J., 119, 1375-1388 (2000). 

12. Hjellming, R. M., et al. Radio observations of the 1985 outburst of RS Ophiuchi. Astrophy. J., 
305, L71-L75 (1986). 

13. Snijders, M. A. J. Ultraviolet Observations of the 1985 Outburst of RS Ophiuchi. in RS Ophi- 
uchi (1985) and the Recurrent Nova Phenomenon, M. F. Bode, Ed., 51-60 (1987). 

14. Mukai, K. & Ishida, M. The Early X-Ray Emission from V382 Velorum (Nova velorum 1999): 
An Internal Shock Model. Astrophys. J., 551, 1024-1030 (2001). 

15. Greiner, J., Orio, M., Schartel, N. XMM-Newton observations of Nova LMC 2000. Astron. 
and Astrophys., 405, 703-710 (2003). 

16. Orio, M. X-Ray Observations of Novae, in IAU Colloquium 194 - Compact Binaries in the 
Galaxy and Beyond, Ed. G. Tovmassian & E. Sion (2004) 



6 



17. Chevalier, R. A., Kirshner, R. P., & Raymond, J. C. The optical emission from a fast shock 
wave with application to supernova remnants. Astrophys. J., 235, 186-195 (1980). 

18. Iijima, T. RS Ophiuchi. IAU Circ, 8675 (2006). 

19. Buil, C. RS Ophiuchi. Cent. Bur. Astron. Tel., 403 (2006). 

20. Fujii, M. Optical spectrum of RS Oph. VSNET-alert, 8869 (2006). 

21. O'Brien, T. J., Lloyd, H. M, & Bode, M. F. An interacting winds model for the X-ray emission 
from V838 HER (nova Herculis 1991). Mon. Not. R. Astron. Soc, 271, 155-160 (1994). 

22. Mason, K. O., Cordova, F. A., Bode, M. F, & Barr, P. Intense Soft X-rays from RS Ophiuchi 
during the 1985 Outburst, in RS Ophiuchi ( 1985) and the Recurrent Nova Phenomenon, M. F. 
Bode, Ed., 167-175 (1987) 

23. O'Brien, T. J., & Kahn, F. D. 1987, Models for the remnants of recurrent novae - II. Dynamical 
effect of radiative heat loss. Mon. Not. R. Astron. Soc, 228, 277-287 (1987). 

24. Taylor, A. R., Davis, R. J., Porcas, R. W., & Bode, M. F. VLBI observations of RS Oph - a 
recurrent nova with non-spherical ejection. Mon. Not. R. Astron. Soc, 237, 81-91 (1989). 

25. Lloyd, H. M., Bode, M. F, O'Brien, T. J., & Kahn, F. D. Bipolar recurrent nova outbursts - I. 
Hydrodynamic models. Mon. Not. R. Astron. Soc, 265, 457-470 (1993). 

26. Hachisu, I., & Kato, M. Recurrent Novae as a Progenitor System f Type la Supernovae. I. 
RS Ophiuchi Subclass: Systems with a Red Giant Companion. Astrophys. J. 558, 323-350 
(2001). 

27. Hamuy, M., et al. An asymptotic-giant-branch star in the progenitor system of a type la super- 
nova. Nature, 424, 651-654(2003). 

28. Wood-Vasey, W. M., Wang, L., & Aldering, G. Photometry of SN 2002ic and Implications for 
the progenitor Mass-Loss History. Astrophys. J., 616, 339-345 (2004). 

29. Wood-Vasey, W. M., & Sokoloski, J. L. Novae as a mechanism for producing cavities around 
the progenitors of SN 2002ic and other SNe la. Astrophys. J., submitted. 

30. Davis, R. J. Radio Observations of RS Ophiuchi. in RS Ophiuchi (1985) and the Recurrent 
Nova Phenomenon, M. F. Bode, Ed., 187-192 (1987). 



Acknowledgements We are grateful for discussions with J. Raymond, C. Rakowski, M. Wood-Vasey, and 
C. Matzner. J.L.S. acknowledges support from the NSF. G.J.M.L. acknowledges support from NASA, plus 
CNPq and FAPESP (Brazil). 

Competing Interests The authors declare that they have no competing financial interests. 



7 



Correspondence Correspondence and requests for materials should be addressed to 
J.L.S. (jsokoloski@cfa.harvard.edu). 



8 



